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Electric  field  data  from  the  double  probe  vector  electric  field  instrument  (VEFI)  on  the  DE  2  spacecraft 
have  been  analyzed  to  determine  the  average  meridional  electric  field  (zonal  ion  flow)  patterns  in  the 
region  between  ±  30  magnetic  latitude  during  solar  maximum  conditions.  Over  300  passes  were  used  to 
compile  the  data  set.  Data  were  projected  to  a  constant  300-km  altitude  and  to  the  equatorial  plane 
assuming  that  the  electric  field  along  the  magnetic  field  was  zero.  The  average  data  set  displayed  a  rapid 
increase  of  the  downward  meridional  electric  field  with  local  time  near  1800  MLT  with  the  higher 
latitudes  seeing  the  change  first  A  secondary  nighttime  maximum  of  this  electric  field  component  was 
observed  post  midnight  with  the  crossover  to  upward  electric  fields  (westward  ion  flow)  occurring 
between  0400  and  0500  MLT.  A  sharp  return  to  near  zero  was  observed  between  1200  and  1300.  Typical 
average  amplitudes  range  between  3  and  6  mV  m.  No  consistent  variations  with  magnetic  activity  were 
observed.  Although  the  daily  variation  in  the  zonal  ion  flow  is  dominated  by  the  diurnal  term,  a  net 
superrotation  is  evident  in  the  harmonic  analysis.  The  superrotation  is  strongest  near  the  equator  and 
decreases  with  latitude,  because  of  the  disturbance  dynamo.  The  higher  order  harmonics  up  through  the 
quatrediurnal  term  are  also  of  significant  magnitude  in  the  analysis  of  the  shape  of  the  daily  variation. 
Close  similarity  is  seen  to  the  zonal  neutral  winds  indicating  that  they  form  the  principal  driving  force. 
The  magnitudes  of  (he  electric  field  derived  ion  drift  are  somewhat  higher  than  the  average  F  region 
neutral  wind  values.  This  along  with  the  higher  order  harmonic  content  argues  for  the  need  to  develop 
fully  coupled  F.  and  F  region  model  depicting  the  ionosphere  and  thermosphere  interactions  in  a 
self-consistent  fashion. 


I.  Introduction 

Observations  of  meridional  electric  fields  or  zonal  plasma 
flow  at  low  and  mid-latitudes  have  in  the  o-st  been  primarily 
made  from  incoherent  scatter  radar  stations  such  as  Jicamarca 
near  Lima.  Peru  [Woodman.  1972;  Fejer  et  al..  1981.  1985], 
Arecibo  in  Puerto  Rico  [Behnke  and  Harper.  1973;  Ganguly  el 
a!.  1987]  and  Saint  Santin  [Blanc  and  Amayene.  1979]  to 
determine  average  ion  drifts  including  seasonal  effects.  Obser¬ 
vations  of  the  drift  of  scintillation  irregularity  structures  has 
also  added  limited  ion  drift  data  at  night  when  spread  F  is 
present  [Abdu  el  al..  1985;  Basil  el  al..  1980].  The  low  latitude 
electric  field  measurements  from  the  DE-2  satellite  [Aggson  ei 
al..  1987]  provide  a  new  more  spatially  extensive  data  source 
with  coverage  at  all  local  times  and  longitudes. 

The  characteristics  of  equatorial  zonal  plasma  flow  have 
been  summarized  by  Fejer  el  al.  [1985]  based  on  Jicamarca 
data.  Westward  drifts  up  to  50  m  s  are  typical  of  daytime 
hours  while  at  night  eastward  drifts  up  to  130  m  s  dominate. 
During  the  day  F  region  dynamo  elfects  [Rishbelh.  1971]  are 
shorted  to  a  large  extent  by  the  high  conductivity  in  the  £ 
regie  i.  hile  after  dusk  the  fields  result  from  the  coupled  £ 
and  F  region  dynamos  [see  Fejer  et  al.  1981],  The  result  is  a 
i.ipij  ineieaxe  in  the  eastwc.d  drift  just  alter  local  sunset.  The 
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equatorial  eastward  flow  increases  sharply  with  increasing 
altitude  below  the  maximum  of  the  F  region,  and  the  altitude 
of  the  steep  gradient  increases  with  increasing  local  time 
[Kiidcki  el  al..  1981;  Aggson  el  at..  1987],  Fejer  et  al.  [1981] 
found  that  the  zonal  drifts  have  considerable  day-to-day 
variability,  but  do  not  show  significant  changes  with  season  or 
magnetic  activity.  This  contrasts  with  the  significant  variations 
with  magnetic  activity  seen  in  the  east-west  electric  field  or 
vertical  drifts  [see  Fejer.  1981],  Fejer  et  al  [1985]  did  find 
increases  in  the  nighttime  maximum  attributable  to  solar  cycle 
variations.  Most  of  the  Jicamarca  data  presented  are  altitude 
averages  from  the  300  to  400  km  altitudes  and  are  treated  as 
height-independent  variations.  Some  latitude  variations  pro¬ 
jected  to  the  equator  as  altitude  variations  may  affect  their 
local  time  averages.  Fejer  et  al.  [1985]  also  found  significant 
fluctuations  on  time  scales  of  5  10  min;  however  they  cau¬ 
tioned  about  over  interpreting  these  variations  due  to  nonop- 
ttmizcd  analysis  procedures. 

The  ion  drifts  (Vi)  can  be  related  to  polarization  electric 
fields  (K)  by 

Vi  =  (E  x  Bl  B1 

Structures  in  the  vertical  profile  of  electric  fields  can  in  turn 
result  from  latitudinal  variations  which  are  projected  along 
equipotential  magnetic  field  lines  to  the  geomagnetic  equator 
Thus  measurements  of  the  meridional  electric  field  obtained 
from  the  polar  orbiting  DE  2  satellite  can  be  projected  to  a 
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latitudinal  profile  at  constant  altitude  or  to  an  altitude  profile 
at  the  equator.  Ayyson  cl  al.  [19871  showed  that  in  the  eve¬ 
ning  the  electric  field  above  the  equatorial  F  region  peak  de¬ 
creased  with  increasing  altitude  or  latitude.  This  effect  had 
been  predicted  by  Anderson  and  Mendillo  (1983)  in  their  cou¬ 
pled  K  and  F  region  model. 

F.  region  dynamo  models  have  been  developed  over  the 
years  to  explain  the  daily  variations  of  ionospheric  currents 
derived  from  ground  magnetometer  measurements.  Despite 
the  apparent  success  of  dynamo  theory  to  model  the  currents, 
there  is  poor  agreement  with  midlatitude  electric  fields  de¬ 
duced  from  ion  drift  measurements  [Richmond  cl  al.,  1976], 
They  added  ion  drag  influences  on  the  wind  to  the  dynamo 
formulations  to  provide  some  improvement.  In  the  more 
recent  models  of  Anderson  aiul  Mendillo  [1983]  and  Takeda 
and  Macda  [1983]  the  self-consistent  coupling  of  the  F  region 
dynamo  [ Rishhclh.  1971 ;  llcclis  cl  al..  1974]  with  the  F.  region 
dynamo  in  which  field  aliened  cui.ent  flow  is  allowed  have 
resulted  in  better  agreement  with  these  data. 

Richmond  ci  al.  [1980]  have  taken  the  radar  data  from 
Arecibo,  Jicamarca.  St.  Santin  and  Millstone  Hill  for  quiet 
periods  Kp  <  2+  and  have  empirically  extrapolated  to  all 
latitudes,  longitudes  and  local  times  using  a  polynomial  ex¬ 
pansion.  Some  seasonal  dependence  is  factored  into  the  model. 
The  amplitude  of  the  zonal  drifts  decreases  with  increasing 
latitude.  A  maximum  of  near  80  m  s  is  seen  near  2100  MLT  at 
10  magnetic  latitude  and  a  mimimum  of  near  40  m  s  west¬ 
ward  drift  is  seen  near  1 1  Of)  MLT.  Zero  crossings  are  at  about 
0500  and  1600  MLT.  The  variation  is  principally  diurnal. 

The  principal  driving  force  for  the  ions  at  low  and  mid 
latitudes  is  the  neutral  wind  Recent  measurements  of  the  F 
region  zonal  neutral  winds  from  DE  2  have  shown  a  basic 
diurnal  variation  with  a  small  superrotation  term  [Wharton  ct 
al..  1984],  Hcrrcro  and  Mayr  [1986]  have  done  a  spectral 
analysis  on  both  the  DE  neutral  wind  data  and  the  Jicamarca 
ion  drift  data  and  found  close  correspondence  on  all  terms 
except  the  superroialion  term.  The  ions  show  larger  superro- 
talion  in  the  equatorial  region,  since  the  westward  daytime 
drift  is  smaller  than  the  eastward  nighttime  drift  due  to  E 
region  loading. 

In  this  paper  we  will  use  the  spatially  more  extensive  DE  2 
electric  field  data  set  to  provide  an  average  picture  of  the 
low -latitude  meridional  electric  fields  or  zonal  plasma  drift. 
These  data  make  it  possible  to  obtain  latitudinal  variations  of 
the  ion  drift.  The  results  will  then  be  compared  to  those  ob¬ 
tained  from  radar  measurements  and  also  to  the  DE  2  neutral 
wind  measurements. 

2.  Data  H  vndi  inc; 

The  Vector  Electric  Field  Instrument  (VEKII  on  DE  2 
[Maynard  cl  al..  1981]  measured  the  electric  field  vector  in  the 
orbit  plane.  For  the  90  inclination  orbit  this  translates  into 
the  meridional  component  in  both  magnetic  and  geographic 
coordinates  at  low  latitudes.  The  DE  2  spacecraft  was  three- 
axis  stabilized  with  the  electric  field  measurement  axes  al 
J_45  to  the  horizontal.  Cylindrical  antennas,  each  11m  m 

I  ’  th.  .  . . ted  m  t  ach  'bi-ectior,  l;,i  each  measurement 

axis.  I  he  sensor  part  of  each  antenna  was  the  outer  2  m 
making  the  baseline  for  the  double  floating  probe  measure¬ 
ment  to  be  2 1 .4  m. 

Previous  studies  with  the  instrument  have  primarily  con¬ 
centrated  on  the  high-latitude  data  where  electric  fields  are 
much  larger  The  measurement  errors  have  been  I  4 mV  m 


due  to  uncertainties  in  constant  offsets  in  the  measured  poten¬ 
tial  difference  usually  referred  to  as  “contact  potential”  [see 
Maynard  ct  al..  1982],  At  high  latitude  the  adjustment  of  a 
constant  offset  for  each  axis  for  the  whole  pass  has  been  made 
assuming  that  the  electric  field  is  near  zero  at  mid-latitude  and 
that  the  integrated  potential  across  the  high-latitude  region 
from  mid-latitude  on  one  side  of  the  globe  to  mid-latitude  on 
the  other  side  is  near  zero.  It  was  puzzling  that  a  different 
“constant"  was  often  needed  for  the  southern  hemisphere  data 
compared  to  the  northern  hemisphere  data  on  the  same  orbit 
This  question  has  been  resolved  by  decomposing  the  low- 
latitude  data  into  components  perpendicular  and  parallel  to 
the  in  orbit  plane  magnetic  field  from  passes  crossing  the  geo¬ 
magnetic  equator.  Since  the  magnetic  field  is  nearly  in  the 
orbit  plane,  one  would  expect  the  parallel  component  of  the 
electric  field  to  be  almost  exactly  zero. 

Figure  I  shows  an  example  of  DE  2  electric  field  data  in 
which  the  contact  potential  adjustment  has  been  determined 
such  that  the  field  aligned  component  E.  (lower  panel)  is  near 
zero  over  most  of  the  regions.  The  constancy  of  a  near  zero  E. 
within  a  broad  latitude  range  and  the  symmetry  of  the  re¬ 
sulting  perpendicular  component  Ex  (middle  panel)  about  the 
geomagnetic  equator  testifies  to  the  correctness  of  the  adjust¬ 
ment. 
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l  ig  1  The  electric  field  measured  b>  DF  2  in  a  magneticalh 
oriented  coordinate  swlem  I  he  bottom  two  panels  show  the  meridi¬ 
onal  component  [hw  (positive  poleward  or  outward)  perpendicular  to 
the  m-plane  magnetic  field  and  the  component  parallel  (/:rl  to  the 
in-plane  magnetic  lield  The  top  panel  shows  the  projection  of  these 
components  along  the  magnetic  field  to  a  constant  altitude  of  MM)  km. 
assuming  that  the  electric  field  along  the  magnetic  lield  is  zero,  versus 
latitude 
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At  each  end  of  the  pass  the  E.  component  deviates  from 
zero.  Examination  of  the  data  in  the  original  measurement 
axes  showed  that  each  change  occurred  in  data  from  one  axis 
alone,  and  that  it  occurred  only  when  that  axis  was  nearly 
aligned  with  the  magnetic  field.  Thus  a  contact  potential 
change  occurs  whenever  an  axis  passes  through  the  magnetic 
wake  of  the  satellite  (i.e..  away  from  the  velocity  vector  and 
along  B).  Each  axis  will  become  nearly  aligned  with  the  mag¬ 
netic  field  once  in  each  low  latitude  passage  explaining  the 
differences  in  contact  potentials  for  northern  and  southern 
high  latitude  regions.  Causes  of  contact  potential  differences 
are  primarily  related  to  asymmetries  in  surface  properties  of 
opposite  sensors  either  from  contamination  or  geometry. 
Thus,  the  above  implies  asymmetric  flow  along  B  and  deposi¬ 
tion  of  contaminant  ions  from  the  spacecraft.  Langmuir 
probes  have  experienced  contaminant  deposition  over  time, 
prompting  the  use  of  heaters  to  clean  the  probe  surface  of 
contaminants  fsee  K rehbiel  el  ai,  1981], 

By  restricting  the  coverage  of  this  study  to  +30  geomag¬ 
netic  it  was  possible  to  automate  the  determination  of  the 
contact  potential  for  each  pass.  The  fact  that  the  V  x  B  elec¬ 
tric  fields  which  have  to  be  subtracted  are  smallest  near  the 
equator  reduces  the  effects  of  attitude  errors.  F  urthermore,  the 
fact  that  the  contact  potential  is  constant  over  this  range 
allows  us  to  reduce  the  measurement  errors  to  +2  rr.V.'m.  Of 
the  over  3(H)  passes  processed  less  than  2° ..  could  not  be 
handled  in  an  automated  fashion  with  a  constant  contact  po¬ 
tential.  These  were  omitted  from  the  study  using  the  criteria 
that  the  \E.\  must  be  less  than  2mV  m  for  the  E„  data  to  be 
valid. 

The  meridional  electric  held  values  measured  at  varying 
altitude  along  the  orbit  track  were  projected  along  the  mag¬ 
netic  Held  to  a  constant  altitude  of  300  km.  This  provides 
assuming  E  •  B  =  0  a  latitude  local  time  map  of  Ex.  Using  the 
same  argument,  the  electric  Held  can  be  projected  to  0  mag¬ 
netic  latitude  giving  the  altitude  profile  of  Ex  at  the  geomag¬ 
netic  equator.  This  is  discussed  in  more  detail  below.  The  top 
panel  in  Figure  I  shows  the  300-km  projection.  The  gap  at  the 
center  of  the  plot  results  from  the  altitude  of  the  equator 
crossing  being  above  300  km.  The  precession  of  perigee 
around  the  orbit  and  the  constancy  of  the  orbit  plane  in  iner¬ 
tial  space  causes  certain  local  times  to  be  only  covered  at 
equinox  and  others  at  solstice,  and  also  restricts  the  data 
coverage  inside  +  15  (latitude  survey)  or  below  375  km  (equa¬ 
torial  projection). 

Figure  2  shows  seven  individual  passes  in  the  0300  0400 
V'  T  slot  projected  to  300  km  constant  altitude.  No  selection 
criteria  were  exercized  except  that  of  local  time.  The  general 
consistency  of  the  data  from  day  to  day  anu  pass  to  pass  is 
evident.  A  montage  of  passes  from  near  2000  MLT  was  pre¬ 
sented  by  .It/c/xan  cl  ai  [1987]  (their  Figure  2).  Again,  the 
overall  consistency  of  the  latitude  profiles  at  a  particular  local 
time  was  demonstrated  from  pass  to  pass.  Typical  profiles 
versus  latitude  are  shown  in  F.gui .-  1  f.-r  evrry  2  hour.”,  rvr 
the  24-hour  day.  Systematic  changes  with  local  time  are  evi¬ 
dent 

For  most  passes  the  electric  field  is  regular  and  well  be¬ 
haved  In  a  few  cases  during  magnetic  storm  periods  spread  F 
irregularities  increased  to  the  point  of  dominating  the  dc  elec¬ 
tric  field  Figuic  4  shows  an  example  near  2000  MLT.  Note 
the  change  in  scale.  Even  though,  the  irregularities  reach 
nearly  20  mV  m.  the  background  electric  field  is  still  in  the 
5  10  mV  m  region  which  is  typical  for  that  local  time.  The 
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Fig.  2  A  collection  of  projected  perpendicular  electric  field  (posi¬ 
tive  poleward  or  outward)  profiles  spanning  about  one  hour  of  local 
lime  illuslraling  the  repeatability  of  the  data  al  a  particular  local  time. 


regions  of  spread  F  are  not  conjugate  in  the  plot  and  should 
not  be  expected  to  be  as  they  are  local  and  not  necessarily 
extended  in  longitude.  Both  time  and  spatial  variations  work 
against  the  probability  of  actually  observing  conjugacy  on  a 
given  pass  The  remainder  of  the  paper  will  concentrate  on  the 
average  electric  field  as  this  type  of  activity  is  not  typical  and 
does  not  grossly  distort  the  averages. 

3.  Average  Patterns 

In  order  to  establish  a  statistical  picture  of  the  meridional 
electric  field  at  low  latitude,  all  projected  data  were  averaged 
in  5  latitudinal  bins  (300-km  projection)  and  25-km  altitude 
bins  (equatorial  projection).  The  results  were  then  combined 
and  averaged  for  each  hour  of  local  time.  Three  hundred  and 
two  passes  '  ere  used  in  this  study  covering  all  local  times. 
Because  of  the  afore  mentioned  nonumtorm  orbital  coverage 
and  a  preponderance  of  passes  in  (he  I9hr  time  frame  used  in 
the  previous  study  [  li/i/son  cl  ai.  19X7]  the  statistics  are  not 
u:;:fe .  ..  ‘TV.b’ •  1  -nows  the  distribution  of  data  points  with 
latitude  north  and  south  of  the  equator.  At  any  given  latitude, 
each  point  represents  a  separate  pass.  Conjugate  symmetry 
has  been  assumed  and  the  northern  and  southern  hemisphere 
data  sets  have  been  combined  for  the  equatorial  projection. 

The  resulting  data  sets  were  averaged  for  each  magnetic 
local  time  latitude  bin  and  are  displayed  in  Figure  s  for  the 
latitudinal  projection  and  Figure  6  for  the  vertical  or  equa¬ 
torial  projection.  The  vectors  are  rotated  into  the  horizontal 
east-west  drift  direction  for  each  of  displaying  and  to  show  the 
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Fig.  3.  Examples  of  latitudinal  profiles  of  the  projected  perpendicular  electric  field  (positive  poleward  or  outward)  for 
every  2  hours  to  illustrate  the  diurnal  variation.  The  variability  of  the  gap  at  the  magnetic  equator  is  a  result  of  the  local 
time  coverage  be»ng  a  function  of  season  and  of  the  change  in  the  latitude  of  perigee  with  time. 


related  motion.  Actual  meridional  electric  field  values  and  the 
standard  deviation  of  each  average  are  given  in  the  appendix. 
Several  characteristics  stand  out  in  the  patterns: 

1.  The  electric  field  increases  rapidly  with  local  time  near 
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Fig.  4.  A  pass  from  a  magnetically  active  period  containing  very 
large  variations  in  the  elcelrie  tieUI  related  .o  spreati  t  The  bolttvm 
panel  shows  the  orbital  \ariation  as  taken  with  time  while  the  top 
panel  shows  the  .a-J  data  projected  to  a  constant  100-km  altitude 
versus  magnetic  latitude  Note  that  even  though  the  variations  are 
large,  the  background  electric  field  is  comparable  to  that  seen  on 
other  passes  at  the  same  local  time. 


1800  hrs  with  the  higher  latitudes  (or  altitudes  in  the  equa¬ 
torial  projections)  seeing  the  increase  first.  This  was  discussed 
by  Aggson  et  al.  [1987], 

2.  The  downward  electric  field  (eastward  drift)  decreases 
near  midnight  and  then  increases  again  before  the  reversal  to 
westward  drift  between  0400  and  0500  MLT. 

3.  Electric  fields  are  upward  corresponding  to  westward 
drifts  in  the  morning  hours. 

4.  A  sharp  return  to  near  zero  electric  fields  is  observed 
between  1200  and  1300  MLT.  The  electric  fields  remain  near 
zero  during  the  remainder  of  the  afternoon. 

5.  Average  electric  field  magnitudes  range  up  to  8  mV  m 
with  3-6  mV  m  being  typical  for  most  local  times. 

6.  The  integrated  areas  under  the  eastward  drift  portions 
of  the  curves  are  larger  than  those  under  the  westward  drift 
portions  indicating  that  a  net  superrotation  term  is  present. 
This  unbalance  is  greatest  at  the  lower  latitudes. 

7.  The  zonal  ion  drift  is  approximately  symmetric  with 
respect  to  the  geomagnetic  equator  (some  seasonal  asym¬ 
metries  are  possible  but  not  separable).  These  features  will  be 
discussed  and  compared  to  other  results  in  subsequent  sec¬ 
tions. 

Several  words  of  caution  relative  to  the  data  set  need  to  be 
emphasized.  The  local  time  coverage  is  determined  by  the 
rotation  of  the  orbit  with  season.  Thus  data  from  1800  2200 
and  0600  1000  MLT  is  primarily  solstice  data,  while  equinox 
data  is  emphasized  in  between  The  solstice  data  is  furfb-r 
biased  toward  December  uala  because  o!  the  abundance  of 
perigee  data  near  the  equator  during  December  1981.  The 
movement  of  perigee  around  the  orbit  plane  results  in  the  lack 
of  low  altitude  data  at  some  local  times  and  an  incomplete 
data  set  with  local  time  at  <  5  magnetic  latitude.  Thus  sea¬ 
sonal  variations  can  not  he  separated  out.  The  high  degree  of 
symmetry  displayed  with  respect  to  the  equator  is  indicative  of 
equipotential  magnetic  field  lines  and  or  little  seasonal  vari¬ 
ation  resulting  from  solstiee.  We  assume  that  the  24-hour  vari¬ 
ations  obtained  here  can  be  used  to  estimate  latitudinal  van- 
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TABLE  I.  Data'Distribution  for  Low-Latitude  Meridional  Electric  Fields 
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ations  of  tidal  components  which  correspond  to  a  mixture  of 
equinox  and  northern  hemisphere  winter  for  the  years  1981- 
1982. 

In  early  analyses  of  the  data  set  it  was  noticed  that  there 


was  a  distinct  variation  with  longitude.  Larger  electric  field 
values  were  found  in  regions  where  the  Earth's  magnetic  field 
was  higher.  All  values  were  therefore  individually  converted  to 
drift  and  reaveraged  to  obtain  the  pattern  shown  in  Figure  7. 


DE-2  Average  Meridional  Electric  Rakts  300  km  Projection 


Fig  5.  The  average  meridional  electric  field  at  300  km  based  on  DE  2  electric  field  measurements  The  vectors  have 
been  rotated  into  the  direction  of  the  corresponding  zonal  ion  drift  for  clarity  of  presentation.  All  vectors  are  therefore 
horizontal  since  the  zonal  component  of  the  electric  field  was  not  measured,  hence  could  not  be  included.  The  scale  in 
electric  field  untts  is  at  the  top  right.  Data  from  all  longitudes  are  average  together  (see  text). 
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DE-2  Average  Meridional  Electric  Field  Equatorial  Projection 


magnetic  local  time 

Fig.  6  The  average  meridional  electric  field  projected  along  the  magnetic  field  to  the  equator  and  shown  versus  altitude 
The  vectors  have  been  rotated  into  the  ion  drift  direction  for  clarity  (see  Figure  5  caption). 


This  pattern  is  invariant  with  longitude  within  the  accuracy  of 
the  data,  affirming  that  the  neutral  wind  driven  ions  move 
consistently  at  all  longitudes  and  that  the  electric  field  longi¬ 
tudinal  variations  are  created  by  longitudinal  changes  in  B  in 
the  expression  E  =  —  V  x  B.  Average  drift  values  are  also 
given  in  the  appendix. 

4.  Comparisons  to  Other  Data 

To  facilitate  comparison  with  other  data  sets,  the  average 
zonal-drift  values  for  the  northern  hemisphere  have  been  plot¬ 
ted  versus  local  time  in  Figure  8.  The  increase  starting  at  1800 
MLT  and  the  secondary  peak  near  0400  MLT  are  clearly 
seen.  Note  that  the  secondary  peak  near  4  hrs  decreases  as  one 
gets  farther  from  the  equator.  The  return  to  zero  magnitude 
near  1200  hrs  is  prevalent  at  all  latitudes. 

Figure  9  shows  the  average  eastward  drifts  measured  by  the 
Jicamarca  backscattcr  radar  facility  [Fejer  et  u/..  1985],  The 
curve  is  for  solar  max  conditions  which  best  corresponds  to 
the  solar  cycle  phase  of  the  DE  data.  The  averaged  zonal  drift 
data  for  10  N  magnetic  latitude  from  DE  is  also  plotted  for 
comparison.  Note  that  the  Jicamarca  data  is  collected  from 
300  to  650  km  in  altitude,  with  350  km  being  the  most  typical 
altitude  for  the  average  values  that  are  plotted.  The  secondary 
peak  near  0400  MLT  is  evident  in  both  data  sets  although  it  is 


more  prominent  in  the  DE  data.  The  Jicamarca  magnitudes 
are  25  to  50", .  less  than  the  DE  values.  The  most  pronounced 
difference  between  the  two  curves  is  the  gradual  return  to 
eastward  drifts  in  the  early  afternoon  with  the  zero  crossing 
near  1600  MLT  observed  at  Jicamarca  compared  to  the  more 
abrupt  transition  to  near  zero  in  the  DE  data.  While  the  range 
of  values  shown  by  Fejer  el  al.  [1985]  for  early  afternoon 
easily  encompasses  zero,  their  average  clearly  remains  west¬ 
ward. 

Ganguly  el  til.  [1987]  have  published  average  zonal  drift 
values  derived  from  backscatter  data  taken  at  Arecibo.  which 
is  al  32  N  geomagnetic  latitude.  Amplitudes  range  from  near 
—  20  to  near  +  30  m/s.  well  below  the  DE  values.  The  Arecibo 
data  is  primarily  from  solar  minimum  conditions.  Fejer  et  al. 
[1985]  have  shown  a  considerable  solar  cycle  effect  in  the 
Jicamarca  data  with  solar  maximum  values  being  50  75"„ 
larger  than  those  from  solar  minimum.  The  secondary  night¬ 
time  peak  near  0400  MLT  is  not  seen  at  Arecibo  consistent 
with  the  DF  data  at  +30  magnetic  latitude,  with  the  transi¬ 
tion  to  westward  drifts  at  Arecibo  occurring  between  1000  and 
0400  MLT  depending  on  season.  The  Arecibo  data  set  does 
show  the  return  to  near  zero  in  the  drift  at  noon  in  the  equin¬ 
ox  data.  The  DE  data  in  the  noontime  and  early  afternoon 
periods  are  equinox  data;  thus  the  two  data  sets  confirm  the 
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DE-2  Average  Zonal  Ion  Flow  300  km  Projection 


MAGNETIC  local  time 

Fig  7  The  average  zonal  ion  flow  from  individual  pass  conversions  of  ihe  meridional  electric  field  lo  zonal  flow.  Data 

from  all  longitudes  are  included. 


noontime  decrease  at  mid-latitudes.  Fejer  et  al.  [1981]  did  not 
find  similar  behavior  with  season  in  the  equatorial  Jimarca 
data  set.  The  Arecibo  summer  data  also  shows  a  crossover  to 
eastward  drifts  near  noon  and  also  a  very  different  early 
morning  behavior.  Their  summer  local  maximum  in  the  west¬ 
ward  drift  near  0300  LT  would  not  be  reflected  in  the  DE 
data,  as  equinox  controls  that  portion  of  local  time. 

Figure  10  compares  the  average  zonal  neutral  wind  vari¬ 
ations  from  the  DE  2  neutral  wind  and  temperature  spec¬ 
trometer  (WATS)  instrument  at  the  geographic  equator  for 
350-km  altitude  with  the  VEFI  data.  The  VEFI  error  bars 
correspond  to  the  statistical  fluctuation  of  the  data.  These 
should  be  added  in  a  RMS  fashion  to  the  previously  men¬ 
tioned  +  2  mV  m  or  ~  +70  m/s  measurement  uncertainty 
after  it  has  been  reduced  by  the  square  root  of  the  number  of 
samples  (see  Table  1).  The  diurnal  pattern  of  this  plot  corre¬ 
sponds  the  closest  to  ion  drift  values  derived  from  the  VEFI 
electric  field  data  (Figure  8).  The  shape  is  very  similar  to  that 
of  the  zonal  ion  flow  with  the  exception  of  the  post  noon 
period.  The  ion  flow  magnitude  is  typically  equal  to  or  larger 
than  the  neutral  wind  magnitude.  Note  that  the  wind  data 
corresponds  the  geographic  equator  while  zonal  ion  drift 
values  arc  keyed  to  the  geomagnetic  equator.  The  close  simi¬ 
larity  of  the  two  data  sets  especially  at  night  again  points  to 
the  neutral  wind  as  a  principal  driving  force;  however,  the 


somewhat  larger  ion  velocity  at  some  local  times  presents  a 
problem  with  this  interpretation  (a  topic  that  will  be  ad¬ 
dressed  further  in  the  discussion). 

5.  Fourier  Analysis 

The  DE  data  are  representative  of  equinox  and  northern 
hemisphere  winter,  and  the  average  24-hour  variations  ob¬ 
tained  show  a  high  degree  of  symmetry  with  respect  to  the 
geomagnetic  equator.  Hence  assuming  that  the  average  24- 
hour  variations  represent  a  mixture  of  these  two  seasons,  an 
estimate  of  the  latitudinal  variation  of  some  of  the  tidal  com¬ 
ponents  of  zonal  ion-drift  can  be  determined  through  Fourier 
analysis.  The  average  daily  variations  of  the  ion  drift  have 
been  Fourier  analyzed  at  each  latitude  interval  from  10  to 
30  in  each  hemisphere.  Figure  1 1  shows  the  average  daily 
variations  obtained  at  15  and  30  magnetic  north  latitudes 
from  fourth-order  Fourier  series  expansion  of  the  correspond¬ 
ing  data  shown  in  Figure  7,  The  tidal  amplitude  spectra  for 
both  latitudes  are  shown  in  the  lower  panel.  Table  2  gives  the 
amplitudes  and  phases  up  through  the  fourth  order  for  each 
latitude  as  well  as  the  neutral  wind  values  based  on  DE  2 
WATS  data  [ Herrera  and  Mayr.  1986]  Superrotation  is  rep¬ 
resented  by  Ihe  prevailing  term  4(0)  in  the  spectra.  Accord¬ 
ingly.  the  amplitude  spectra  show  that  4(0)  at  15  is  signifi¬ 
cantly  larger  than  4(0)  at  30  showing  that  the  ions  superro- 
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6.  Discussion 

A  definite  ion  superrotation  is  an  obvious  feature  of  the 
data.  An  examination  of  Figure  8  dearly  shows  more  integrat¬ 
ed  area  under  the  curves  above  zero  than  below  for  all  except 
the  highest  latitude  curves.  The  harmonic  analysis  zeroth 
order  term  displayed  in  Figure  12  indicated  an  average  value 
of  48  m  s  for  the  region  inside  +  15  .  Wharton  et  al.  [1983] 
found  a  neutral  superrotation  of  10  m  s  from  the  DE-2  WATS 
data,  while  Woodman  [1972]  detected  an  ion  superrotation  of 
about  30  m  s  in  the  Jicamarca  data.  Herrera  and  Mayr  [1986] 
showed  that  the  ion  superrotation  amplitude  from  Jicamarca 
lata  increased  with  solar  activity,  with  33  m  s  being  appli¬ 
cable  to  the  1978  1981  period.  The  DE  data  corresponds  to 
higher  solar  activity,  being  more  nearly  centered  on  solar 
maximum 

The  latitudinal  variation  of  the  ion  superrotation  is  perhaps 
the  most  significant  feature  of  these  data.  Examination  of  the 
individual  plots  for  each  magnetic  latitude  in  Figure  8  (and 
also  in  Figure  12)  shows  that  the  ion  drifts  are  strongly  east¬ 
ward  at  low  latitudes  in  the  postmidnight  to  morning  sector, 
and  at  latitudes  near  30  the  corresponding  drift  is  small  and 
possibly  westward.  It  is  known  that  the  ionosphere  subrotates 
during  strong  geomagnetic  disturbances,  as  noted  and  ex¬ 
plained  by  Rhine  and  Richmond  [1980]  in  terms  of  the  distur¬ 
bance  dynamo.  In  the  absence  of  magnetic  disturbances  there 
still  exists  a  quiescent  level  of  auroral  heating.  Indeed,  such  a 
heat  source  has  been  regularly  hypothesized  in  order  to  under¬ 
stand  thermospheric  circulation  phenomena  [Dickinson  et  al.. 
Iv75.  1977],  A  constant  auroral  heat  source  would  produce 
winds  initially  moving  equatorward.  but  turning  westward  due 
to  the  strong  Coriolis  force  component  at  high  latitudes.  Pre¬ 
sumably  this  westward  flow  develops  in  such  a  way  that  it 
decreases  with  latitudinal  displacement  from  the  auroral 


i 


region.  Thus  when  interacting  with  the  diurnal  and  prevailing 
zonal  drifts  excited  at  the  equator,  the  westward  flow  originat¬ 
ing  in  the  polar  region  would  impose  a  tendency  to  buck  the 
net  eastward  flow  in  the  ions,  thus  producing  the  gradual 
decrease  in  superrotation  with  increasing  latitude. 

A  striking  feature  of  the  daily  variation  in  Figure  8  is  the 
abrupt  abatemenl  of  the  westward  flow  that  occurs  near  local 
noon.  A  similar  abatement,  though  not  so  abrupt  occurs  at 
equinox  and  in  summer  in  the  Arecibo  data  [Gaiu/uly  el  al.. 
1987]  an  .  also  in  the  St.  Santin  observations  [ Blanc  and 
I nteyenc.  1979],  thus  establishing  this  feature  as  characteristic 
of  mid-latitudes  at  equinox.  This  is  different  from  the  equa¬ 
torial  Jicamarca  data  which  has  a  more  gradual  transition  (sec 
Figure  9).  Figure  8  shows  that  this  feature  is  present  at  all 
latitudes  in  the  VEFI  data  from  10  to  30  .  As  discussed  in 
section  3.  the  VEFI  data  around  local  noon  where  the  abate¬ 
menl  is  seen  are  equinox  data  corresponding  to  the  Arecibo 
results.  Since  no  longitude  variations  arc  detectable  and  since 
the  sample  distribution  is  relatively  uniform,  the  VEFI  data 
suggests  that  the  Ganguly  el  al  equinox  feature  extends 
toward  the  equator.  This  extension  into  the  region  where  the 
DE  2  WATS  F  region  nctural  wind  data  docs  not  sec  the 
abatemenl  suggests  that  it  may  be  E  region  controlled. 

The  harmonic  analyses  (see  Figure  1 1  and  Table  2)  indicate 
that  the  diurnal,  ter-diurnal  and  quatre-diurnal  terms  all 
maintain  relatively  constant  magnitude  (peaking  slightly  near 
the  equator)  and  phase  with  latitude,  while  the  semidiurnal 
term  decreases  as  the  equator  is  approached.  Herrera  and 
Mayr  [1986]  also  found  a  smaller  semidiurnal  term  compared 
to  the  ter-diurnal  term  in  their  harmonic  analysis  of  the  DE  2 
WATS  neutral  wind  data  for  the  geographic  equator.  The 
shifting  of  the  phase  of  the  larger  semidiurnal  term  to  earlier 
local  times  at  higher  latitudes,  its  reduced  magnitude  al  lower 
latitudes,  and  an  increase  in  phase  with  latitude  for  the  higher 
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Tidal  Amplitudes  and  Phases  at  15*  and  30*  MLAT 


Fig.  II  The  harmonic  analysis  of  the  zonal  ion  flow.  The  bottom  panel  represents  the  amplitude  versus  harmonic 
number  with  I  referring  to  the  34-hour  diurnal  variation  The  top  panel  displays  the  resulting  curves  for  13  and  30  north 
latitude  data  when  the  terms  through  fourth  order  are  included. 
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order  harmonics  apparently  results  in  the  decrease  in  the  sec¬ 
ondary  nighttime  peak  near  0400  MLT  with  increasing  lati¬ 
tude  (compare  the  fits  for  15  and  30  in  the  upper  panel  in 
Figure  II.  Table  2.  and  the  corresponding  data  curves  in 
Figure  8).  The  origin  of  the  higher-order  harmonics  can  be  in 
the  form  of  the  basic  solar  forcing  function  or  from  modal 
coupling  due  to  nonlinearities  in  the  system  and  to  interhe¬ 
misphere  coupling.  Herrera  el  al.  [1983]  found  significant  am¬ 
plitudes  in  the  higher  order  harmonics  in  their  tidal  analysis  of 
the  thermospheric  midnight  temperature  maximum  which  is 


related  to  the  driving  term  for  these  drifts  through  the  neutral 
density  variations.  Again  in  that  data  set  the  ter-diurnal  am¬ 
plitudes  were  sometimes  larger  than  the  semidiurnal  ampli¬ 
tudes.  Caution  in  over-interpreting  the  higher  order  harmon¬ 
ics  should  be  taken  in  light  of  their  low  magnitude  and  the 
seasonal  variability  in  the  Arecibo  data  [Ganguly  el  al,  1987] 
which  is  obviously  mixed  into  in  this  analysis 

larhes  anj  Gillette  [1982]  provided  a  compendium  of  theo¬ 
retical  atmospheric  tidal  structures  based  on  diurnal  and  semi¬ 
diurnal  terms  The  upper  thermosphere  diurnal  components 
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TABLE  2.  Harmonic  Analyses  Amplitudes  and  Phases 


Magnetic  Latitude, 
deg 

A(0), 

m/s 

All) 

(PHI), 
m/s  (hour) 

A(2) 
(PH2), 
m/s  (hour) 

AO) 

(PH3), 
m/s  (hour) 

A(4) 
(PH4), 
m/s  (hour) 

Zona 1  ton  Drift  Velocity 

-30 

19 

112  (20.8) 

32  (0.3) 

38  (16.3) 

33  (13.2) 

-25 

32 

116  (21.2) 

37  (-0.3) 

37  (16.8) 

30  (14.0) 

-20 

47 

120  (21.7) 

25  (-0.1) 

37  (15.4) 

38  (13.4) 

-15 

54 

128  (21.9) 

12  (1.3) 

41  (14.5) 

40  (12.6) 

-10 

59 

134  (22.3) 

22  (7.4) 

46  (13.4) 

50(11.8) 

10 

46 

141  (22.3) 

16  (5.4) 

41  (14.4) 

41  (12.8) 

15 

32 

133  (21.9) 

12(1.6) 

36  (15.3) 

36  (13.4) 

20 

17 

122  (21.5) 

23  (0) 

29  (17.4) 

28  (12.8) 

25 

4 

117(21.2) 

28  (-1.0) 

33  (18.1) 

26  (12.4) 

30 

7 

120  (20.8) 

34  (-1.7) 

35  (18.4) 

26  (12.6) 

Wind  Velocity  [From  Herrera  and  Mayr% 

1986] 

Geographic 

10 

119(22.7) 

18  (15.4) 

32  1 1 1.0) 

20  (14  0) 

Equator 

are  of  lower  magnitude  (nearer  70  m/s)  but  within  about  an 
hour  in  phase  depending  on  the  altitude  (the  theoretical  phase 
is  equal  to  our  average  value  at  180  km  and  about  one  hour 
less  at  300  km).  Their  solar  semidiurnal  tide  is  of  the  order  of 
30  m,s  with  nearly  the  same  phase  variation  that  is  observed 
h«*-e  with  latitude  for  300-km  altitude.  However,  the  ampli¬ 
tude  variation  with  latitude  of  the  F  region  semidiurnal  term 
was  either  constant  or  peaking  near  the  equator.  The  E  region 
semidiurnal  term  is  larger  at  mid-latitudes  with  an  8-10  hour 
phase  difference.  In  general,  the  Forbes  and  Gillette  tidal 
structure  is  diurnally  dominated  in  the  upper  thermosphere 
with  that  domination  increasing  with  latitude.  Note  that  the 
semidiurnal  term  is  the  principal  term  in  the  E  region  (see  also 
Richmond  el  al.  [1976]).  The  diurnal  domination  of  the  VEFI 
derived  ion  drifts  with  similar  phase  to  the  theoretical  struc¬ 
ture  indicates  that  the  F  region  dynamo  plays  a  major  role  in 
the  electric  field  structure.  The  variation  of  the  semidiurnal 
term  and  the  coupling  of  the  higher  harmonics  is  indicative 
that  the  coupling  to  the  E  region  may  also  be  important 
factor,  especially  near  noon  were  a  strong  mid-latitude  E 
region  semidiurnal  term  could  play  a  large  role  in  decoupling 
the  ion  drifts  from  the  F  region  neutral  winds. 

The  largest  difference  between  the  electric  field  derived  drift 
velocities  presented  here  and  previous  results,  both  theoretical 
and  experimental,  is  in  the  amplitude.  The  generally  higher 


DE-2  Zonal  Ion  Flow 
Harmonic  Analytic 
Prevailing  (Super  Rotation)  Term 


Fir  12  The  variation  of  ihe  prevailing  or  super  rotating  zonal  ion 
flow  with  magnetic  latitude. 


values  reported  here  compared  to  the  radar  measurements  can 
to  some  degree  be  explained  by  differences  in  solar  cycle,  the 
day-to-day  variability  in  the  data,  and  limitations  in  the  sizes 
and  overlap  of  the  data  sets.  The  hardest  difference  to  explain 
is  the  comparison  with  the  WATS  data  in  Figure  10.  Both 
data  sets  are  from  the  same  satellite,  although  Ihey  were  hand¬ 
led  differently  (geomagnetic  focus  versus  geographic  focus)  and 
the  pass  selection  was  different.  As  can  be  seen  from  figure  10. 
the  largest  discrepancy,  irrespective  at  the  error  bars,  between 
VEFI  and  WATS  data  set  is  that  the  night  time  averaged 
plasma  drifts  from  this  study  up  to  100  m/s  larger  than  the 
neutral  wind  velocities  as  averaged  in  the  WATS  study 
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l  ig.  13  The  diurnal  zonal  ion  flow  and  its  phase  versus  magnetic 
latitude. 
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[  Wharton  et  al.,  1982].  The  largest  differences  occur  2  hours 
after  sunset  and  2  hours  before  sunrise.  However,  the  periods 
of  largest  discrepancy  are  also  periods  at  the  largest  error 
bars.  (±100  m/'s).  More  exact  comparison  between  the 
averages  of  this  study  and  that  of  Wharton  et  al.  should  prop¬ 
erly  include  in  the  future,  orbit  by  orbit  comparisons  rather 
than  simple  averages.  In  the  Fourier  domain,  as  can  be  seen 
from  Table  2.  the  primary  difference  show  up  in  a  larger  su- 
perrotation  term  and  a  slightly  larger  diurnal  term.  We  would 
have  anticipated  that  the  nighttime  ion  drifts  would  have  been 
comparable  to  the  winds  (not  larger),  while  the  daytime  drifts 
would  have  been  smaller  than  the  winds  (not  comparable), 
based  on  the  daytime  shorting  of  the  F  region  dynamo  [Rish- 
ht’lh.  1971]  by  the  daytime  F.  region  conductivity.  This  mecha¬ 
nism  creates  a  larger  superrotation  term  but  requires  a  smaller 
diurnal  term  While  it  was  noted  by  Fejer  el  al.  [1985]  that 
ground-based  neutral  wind  measurements  were  typically  less 
than  the  Jicamarca  ion  drift  measurements,  they  also  pointed 
out  as  the  cause  'hat  the  ground  based  wind  measurements 
were  related  to  a  lower  altitude  (near  180  km).  A  similar  mag¬ 
nitude  problem  exists  in  comparison  to  the  theoretical  calcula¬ 
tion  of  Forhe s  mul  Gillette  [I982|  in  that  their  magnitude  for 
the  diurnal  term  is  less  The  inaccuracies  in  the  measurement 
(  +  2  mV  m  or  about  +  70  m  s)  when  added  in  RMS  fashion  to 
the  statistical  variations  shown  in  Figure  10  allow  most  of  the 


WATS  curve  to  lie  within  the  error  bar  envelope  or  at  larger 
magnitudes;  however,  the  basic  trend  for  higher  drift  values  is 
evident.  The  validity  of  this  trend  can  only  be  examined 
through  self  consistently  coupled  calculations  of  the  iono¬ 
sphere  and  the  thermosphere  involving  both  the  E  region  and 
the  F  region,  since  sources  and  sinks  in  both  regions  are  in¬ 
volved. 

7.  Summary 

Electric  field  data  from  over  300  equatorial  passes  of  the 
DE  2  satellite  have  been  analyzed  to  produce  average  meridio¬ 
nal  electric  held  or  zonal  ion  How  patterns  for  the  region 
between  +30  geomagnetic  latitude.  Several  characteristics 
stand  out : 

1.  The  downward  meridional  electric  field  increases  rapily 
with  local  time  near  1800  MLT  with  the  higher  latitudes  (or 
altitudes  in  the  equatorial  projection)  seeing  the  change  first 
(see  also  Ayi/son  et  al.  [1987]). 

2.  The  downward  electric  field  (eastward  drift)  decreases 
near  midnight  and  then  has  a  secondary  maximum  before 
reversing  to  upward  electric  field  (westward  drift)  between 
0400  and  0500  MLT. 

3.  A  sharp  return  to  zero  electric  field  is  observed  between 
1200  and  1300  MLT  and  the  electric  field  remains  near  zero 
throughout  the  remainder  of  the  afternoon 

4.  Average  electric  field  magnitudes  range  up  to  8  mV  m 
with  3-6  mV  m  being  most  typical. 

5.  In  general,  the  small-scale  spatial  variations  seen  in  in¬ 
dividual  passes  are  smaller  in  magnitude  than  the  large-scale 
electric  fields  except  during  periods  of  large  magnetic  activity. 

From  harmonic  analysis  of  the  daily  variation  at  each  lati¬ 
tude.  the  following  is  inferred: 

1.  A  net  superrotation  exists  which  is  strongest  at  the 
lowest  latitudes 

2.  The  basic  variation  of  the  ion  flow  is  dominated  by  the 
diurnal  term  which  is  nearly  constant  with  latitude. 

3.  The  semi-diurnal  term  is  weakest  at  the  equator  indicat¬ 
ing  a  mixing  of  t  and  F  region  sources. 

4.  The  third-  and  fourth-order  tidal  terms  are  also  signifi¬ 
cant  and  apparently,  along  with  the  larger  mid-latitude  E 
region  semidiurnal  term,  are  related  to  the  rapid  return  to 
zero  post  noon.  This  feature  in  the  data  is  also  seen  in  the 
Arecibo  and  St.  Santin  drift  data. 

No  significant  variation  in  the  meridional  electric  fields  with 
magnetic  activity  was  observed  consistent  with  previous  re¬ 
sults.  Seasonal  dependencies  are  inferred  but  not  separable 
with  this  data  set.  An  apparent  longitudinal  dependence  in  the 
electric  field  was  shown  to  have  been  the  result  of  the  interac¬ 
tion  of  a  longitudinal-dependence-free  ion  velocity  and  the 
Earth's  magnetic  field. 

The  amplitudes  of  the  electric  field  are  somewhat  higher 
than  would  be  expected  when  compared  to  the  neutral  wind 
measurements  from  l)F.  2  based  on  present  understanding  of 
the  /.  and  F  region  dynamos  The  general  pattern  consistency 
with  both  neutral  wind  measurements  and  with  theoretical 
analyses  tends  to  confirm  the  winds  as  the  source  of  the  elec¬ 
tric  fields  There  is.  however,  no  self  consistent  calculation  that 
couples  the  I.  and  /  region  sources.  Such  a  calculation  will  be 
necessary  to  understand  the  apparent  amplitude  discrepancy 
and  to  establish  the  means  by  which  the  higher  harmonics  are 
created  and  the  role  that  they  play. 

Am  ndix 

The  values  for  the  average  patterns  presented  in  F  igures  5  -7 
are  given  in  Tables  AI-A4  Table  Al  shows  the  meridional 
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electric  fields  perpendicular  to  the  magnetic  field  versus  mag¬ 
netic  latitude  for  the  300-km  projection.  Tables  A2  and  A3 
display  the  vertical  electric  field  versus  altitude  at  the  mag¬ 
netic  equator  for  the  projection  to  the  equatorial  plane.  Table 
A4  lists  the  average  zonal  ion  drift  velocities  derived  from  the 
electric  field  data  versus  latitude  for  the  300-km  projection. 

Each  table  also  presents  values  for  the  standard  deviation  of 
the  data  set  for  each  point.  The  sample  size  at  each  point  in 
the  300-km  projection  is  given  in  Table  1.  Since  the  data  set  is 
so  small  at  the  equator  because  of  orbit  coverage  (see  main 
text),  no  equatorial  samples  are  included.  The  coverage  for  the 
equatorial  projection  is  similar,  with  the  higher  altitudes  cor¬ 
responding  to  the  higher  latitudes  and  the  lower  altitudes  cor¬ 
responding  to  the  lower  latitudes  Where  the  standard  devi¬ 
ation  is  omitted  in  the  tables  represents  data  sets  of  only  one 
value. 

Acknowledgments.  The  Editor  thanks  J.  H.  A.  Sobral  and  another 
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